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ABSTRACT 

We study the growth of dark matter halos in the concordance ACDM cosmology using several N-body 
simulations of large cosmological volumes. We construct merger trees from the Millennium and Millennium-II 
simulations, covering the ranges 10"^- 10'^ Mq in halo mass and 1 - 10^ in merger mass ratio. Our algorithm 
takes special care of halo fragmentation and ensures that the mass contribution of each merger to halo growth 
is only counted once. This way the integrated merger rate converges and we can consistently determine the 
contribution of mergers of different mass ratios to halo growth. We find that all resolved mergers, up to mass 
ratios of 10^ : 1, contribute only « 60% of the total halo mass growth, while major mergers are subdominant, 
e.g. mergers with mass ratios smaller than 3 : 1 (10 : 1) contribute only « 20% (« 30%). This is verified with 
an analysis of two additional simulation boxes, where we follow all particles individually throughout cosmic 
time. Our results are also robust against using several halo definitions. Under the assumption that the power- 
law behaviour of the merger rate at large mass ratios can be extrapolated to arbitrarily large mass ratios, it is 
found that, independent of halo mass, w 40% of the mass in halos comes from genuinely smooth accretion 
of dark matter that was never bound in smaller halos. We discuss possible implications of our findings for 
galaxy formation. One implication, assuming as is standard that the pristine intergalactic medium is heated and 
photoionized by UV photons, is that all halos accrete > 40% of their baryons in smooth "cold" T > 10* K gas, 
rather than as warm, enriched or clumpy gas or as stars. 

Subject headings: cosmology; theory — dark matter — large-scale structure of universe — galaxies: evolution 
— galaxies: formation 



1. INTRODUCTION 

The way galaxies gain their mass affects almost ev- 
ery aspect of galaxy evolution. The distinction between, 
e.g., accretion of gas versus stars, spherical versus filamen- 
tary accretion, or clumpy versus smooth accretion will re- 
sult in very different star-formation histories, colours, mor- 
phologies, angular momentum contents and sizes. Merg- 
ers are believed to play an important role in the evolu- 
tion of galaxies, in particular of elliptical galaxies via the 
morphological transformation fro m disk-dominated galaxies 
to spheroids (e.g. [To omre 1977; Barnes & Hernquist 1992; 
Naab & Bur ker^ 12^)3; Naab et al. .2006: .Cox et al 2006: 
Hopkins et al. 2008; Naab et al. 2001 IConselice et 10120031: 
Bell et al.ll2006[iLotz et alji200 8). Both theoretical and obser- 
vational work have also emphasised the importance of smooth 
accretion of gas, in particular for the buildup of massive 
galaxies at high redshift and for the subs equent evolution of 
disk galaxies (e.g . [White & Frenk 199ll ; iMuraliet alj|2002|; 
Kere§ et all '2005'; iQcvirk et al. 2008; 'Dekel et al.' '2009b'; 
Goerdt et a l. 2009; Daddi et al. 2007; Forster Schreiber et al. 
200^,ConseHce & Arnold 2009; Kauff mann et al.l201Qi) . The 
gas can be accreted in a spherically symmetric mode of cool- 
ing halo gas or in a filamentary mode directly from the cosmic 
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web. 

In the cold dark matt er paradigm for structure formation 
(iBlumenthal et al.lll984l) . galaxy forma tion is closely related 
to the formation of dark matter halos d White & Reesl[l978.) . 
albeit in a complex way. Baryons follow the flow of the 
gravitationally-dominating dark matter and fall into dark mat- 
ter halos, and galaxy mergers follow the mergers of their host 
dark matter halos. Understanding the process of halo mass 
assembly is an important step towards a better understanding 
of galaxy formation. In this paper we use high resolution dark 
matter simulations to study the relative importance of mergers 
versus smooth accretion for the buildup of halos and galaxies. 

The non-linear nature of the evolution of gravitational in- 
stabilities into virialised structures makes N-body simulations 
the most reliable tool available for studying the mass buildup 
of dark matter halos. Massive dark matter particles in N- 
body simulations are believed to give a good representation 
of the coarse-grained phase-space structure of real dark mat- 
ter. However, the interpretation of dark matter simulations is 
subject to some freedom and uncertainty in subsequent steps 
in making the connection to galaxies, such a s the definition 
of a dark matter ha lo (e.g. White 2001, 2002; ICohn & White! 
2008; nnker et alj|2008n and algorithms of merger tree con- 
struction. 

One of the largest simulations used so far for studying the 
growth of dark matter halos is the Millen nium Sirnulation 
ilSpringeletai] 12005: hereafter MS). FakhoiM&Mi jIOOSl) 
analysed the MS and found that the dark matter halo merger 
rate has a nearly universal form that can be separated into 
its dependencies on mass ratio, descendant mass, and red- 
shift. They presented three algorithms for merger tree con- 
struction th a t res ult in merger rates differing by « 25%. In 
iGenel et al] (1200 9) (hereafter G09) we introduced a novel 
mer ger tree construction al gori thm we termed "sp litting" (see 
also lFakhouri & Mall2009l and lMalleret"ani2006l) that incor- 
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porates the complicated process of halo fragmentation and re- 
merging such that a given pair of halos is never considered to 
merge more than once. This resulted in a different set of pa- 
rameters for the Fakhouri & Ma (2008) global fitting formula 
for the merger rate extracted from our trees of the MS. 

Thanks to ever increasing computation power, the dynamic 
range of N-body simulations increases too. As a result, sev- 
eral recent work s have studied the gr owth rate of dark mat- 
ter halos. When IStewart et al.l (l2008h integrated the merger 
contribution and extrapolated to include all unresolved merg- 
ers, they found that only « 50 - 70% of the final mass of 
halos was asse mbled by mergers. Rather dissimilar results 
were found by iMadau et al.l (l2008h . who investigated the for- 
mation ofonejialojna cosmological 'zoom-in' simulation, 
and bv lAngulo & Whita (1201 0!) . who investigated the growth 
of dark matter h alos using high-resolution extended Press- 
Schechter (E PS; iPress & SchechteJITOTl iBond et all fT99Tt 
lBoweij[l99l!) trees. These works concluded that the mass ac- 
cretion of dark matter halos is largely dominated by merg- 
ers. In this paper we perform a similar analysis to that of 
I Stewart etal. (2008 ) with a larger dynamic range using a com- 
bination of the MS and the higher reso lution Millennium-II 
Simulation (iBoylan-Kolchin et al.ll2009l: hereafter MS2), and 
emphasise the importance of choosing an appropriate merger 
tree construction algorithm. We show that all resolved merg- 
ers contribute < 60% to halo mass growth, and suggest that 
« 40% of the accretion rate may be genuinely smooth. We 
also provide further support for these results by directly fol- 
lowing the histories of dark matter particles in two smaller 
cosmological boxes, and make detailed comparisons with the 
previous results mentioned above. 

This paper is organised as follows. In ^2.1l we review the 
Millennium simulatio ns, str ucture identification and merger- 
tree construction. In ^2.2.1l we provide a fittin g function that 
describes the halo merger rate, and in ^2.2.21 we discuss the 
contribution of mergers and smooth accretion to halo mass 
bui ldup . In ^2.2.3l we compare our results to the EPS model. 
In 33.1l we describe our direct analysis of dark matter particles 
histories in two additional simulations, and in 33.2l we present 
the results of that analysis. We discuss implications of our 
results to galaxy formation in ^ and compare our results to 
previous work in ^ We summarise in ^ 

2. ANALYSIS OF MERGER TREES 

2.1. Method 

The MS is a cosmological N-body simulation that follows 
2160^ dark matter particles, each of mass 8.6 x 10^/i~' M0, 
in a periodic box of SOO/i"' Mpc on a side. The cosmology 
is ACDM with n^ = 0.25, Qa = 0.75, ilb = 0.045, h = 0.73, 
n = 1 and erg = 0.9. The MS2 uses the same cosmology and 
follows the same number of particles, but in a box 5 times 
smaller on a side. Thus, the MS2 particle mass is 125 times 
smaller, i.e. 6.885 x lO^/z"' Mq. 

Structure identification in the simulations proceeds in 
two steps. First, the Friends-Of- Friends (EOF) algorithm 
dPavis et al]|1985l: with a linking length parameter b = 0.2) 
creates at every snapshot a catalogue of EOF groups that are 
considered to repres ent dark matter halo s. Second, the al- 
gorithm SUBEIND fSpringel etal .11200 Ih identifies subhalos 
inside EOF groups by finding gravitationally self-bound col- 
lections of particles around maxima in the smoothed density 
field. The terminology is such that even smooth EOF groups 
with no identified substructures contain one subhalo (often 
referred to in the literature as the "main" or "background" 



subhalo) when they are self-bound. Therefore, EOF groups 
that contain zero subhalos are not gravitationally bound and 
are subsequently dropped from the analysis. Publicly avail- 
able subhalo merger trees^ were constructed by finding a sin- 
gle descendant for each subhalo, a procedure in which the 
EOF groups themselves played no role. However, if the halo 
merger rate is to be studied, then different types of merger 
trees need to be built, in which each node is a halo rather 
than a subhalo. Such a well-defined halo merger tree can be 
constructed if each EOF group is given exactly one descen- 
dant. In practice, however, choosing the correct descendant is 
not trivial because EOF groups not only merge, but may also 
fragment back into several groups. 

In most fragmentation events, the subhalo (or group of sub- 
halos) that left its original EOF group becomes a new distinct 
EOF group. Regardless of whether the two merge back to- 
gether (as is usually the case) or not, a sequence of a merger 
followed by a fragmentation introduces artificial effects to the 
statistics of the merger rate and mass accretion. First, such 
a sequence of events is recorded as a merger even though 
the two progenitors end up as two distinct halos at the end 
of the merger-fragmentation sequence. Thus, the merger rate 
is overestimated. Second, a mass accretion rate of M small / ^t 
is attributed to the false merger (Msmaii being the mass of the 
smaller companion and Af the time difference between snap- 
shots), thus the contribution of mergers to halo growth is over- 
estimated too. Third, if no "fragmentation rate" is quantified 
in parallel with the merger rate, all mass changes not associ- 
ated with mergers are considered "smooth", and so a negative 
contribution of -Msmaii/^t is added to the smooth accretion 
component at the time of the fragmentation. Thus, the contri- 
bution of smooth accretion to halo growth is underestimated. 
Additionally, in trees where each halo is allowed to have one 
descendent at most (the standard case in the literature), the 
smaller product of the fragmentation has no progenitor, and 
so its main progenitor track is 'snipped' and all the informa- 
tion of its past formation history is in practice erased. 

Therefore, an algorithm that builds fragmentation-free 
merger trees is needed to quantify correctly both the merger 
rate and the relative contributions of mergers and smooth ac- 
cretion to halo mass growth. In G09 we presented such an 
algorithm and built new trees for the MS. Here we implement 
this algorithm on the MS2 as well. We construct the trees by 
splitting certain EOF groups, those that will suffer a fragmen- 
tation in the future, into several fragments. All the new frag- 
ments that our algorithm creates, as well as untouched EOF 
groups, are considered hereafter simply as "halos". A unique 
descendant is found for every halo, so that the merger tree is 
well defined. More details and motivation for our "splitting" 
algorithm and comparisons to other algorithms ("snipping", 
"stitching" and variants, as well as c ombinations, thereof) are 
presented in G09 and in Fakhouri & Ma (2009). It is worth 
noting that already Mailer et al.. (2006) used a combination of 
"stitching" and "splitting" methods for their N-body/SPH sim- 
ulation to obtain a fragmentation-free merger tree of galaxies. 
The halo merger trees we built (for the MS) are available at 
http://www.mpe.mpg.de/ir/MillenniumMergerTrees/ 

We derive the merger rate per descendant halo per mass ra- 
tio X per unit time to w 1. 69 /D(z), which is the natural time 
variable in the EPS model. Here, D{z) is the linear growth 
rate of density fluctuations, and w is estimated using the 

^ Structure catalogues and merger trees were made pub lic by the Virgo 
Consortium: |http://www.mpa-garching.mpg.de/millenniuin] 
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Fig. 1. — The merger rate per unit time from tlie MS (triangles) and the MS2 (circles). The different masses and redshifts are indicated in each panel. Our 
global fit (Equation[l}, which has an asymptotic power-law index of b = —0.3 at high mass ratios, is shown by red solid lines. In each panel constant power-law 
indices of {—0.4,— 0.3,— 0.2,— 0.1} (from top to bottom) ar'e shown by dashed black lines, and demonstrate that b < —0.2 is never a good fit to the data at z > 0.5. 
At redshifts z < 0.5 our splitting algorithm cannot work properly due to the proximity to the end of the simulation: spurious mergers exist in the trees and the 
power-law index h increases towards zero. Therefore results for low redshift are not shown here, and were not considered for the global fit. 



iNeistein & Dekell (l2008h approximation. In our bookkeeping 
a merger between two halos of masses Mi and M2 < Mi is 
recorded as a merger at mass Mi +M2 with ratio x = M1/M2. 
We define the halo mass as the mass of all particles gravita- 
tionally b ound to it, i.e. the su m of its subhalo masses (see 
G09 and ,Fakhouri & MallMoh . 

2.2. Results from the merger trees 

2.2.1. The merger rate 

We find that the merger rates in the MS and MS2 agree very 
well in the range of overlap (M w 10^^- W^Mq). While 
the MS2, simulating a smaller volume, has worse statistics in 
that range, combining it with the MS provides a much larger 
dynamic range (see below). We fit the merger rate using the 
fitting formula introduced by Fakhouri & Ma (2008) (albeit 
with our mass ratio variable x = I/O- The new best-fitting 
parameters we find are only slightly different from those we 
found for the MS alone in G09. The fitting formula is 



1 



dNr 



merger 



A^desc-lialo dujdx 



(x,z,M)=AM'^2x''exp((x/xr) (1) 



where M12 = M/IO'^ M0. Our best-fitting parameters for the 
combination of the two simulations are; A = 0.065, a = 0.15, 
/7 = -0.3,i = 2.5and7 = 0.5l 



* With the mass ratio definition used bv lFakhouri & M j J2008I) our param- 
eters con'espond, in their notation, to: A = 0.065, « = 0.15, /3 = — 6— 2 = —1.7, 
i=l/x = 0.4, 7 = 0.5, 77 = 1 and M = 10'- Mq . 



Figure [T] shows a few examples of the merger rate for dif- 
ferent halo masses and redshifts and the corresponding fits. 
We find our fit to be appropriate for the whole mass range 
probed by the simulations at redshifts 0.5 ;$ z ;$ 5. At z > 5 
the redshift and mass dependencies become stronger, and we 
do not attempt to fit that regime. At z < 0.5 the fit fails due 
to the proximity to the end of the simulations. That is, the 
subhalo disruption time at z < 0.5 becomes comparable to the 
lookback time corresponding to that redshift, thus we cannot 
identify all the artificially connected halos that would have 
fragmented had the simulations run past z = 0. Figure[T]shows 
that the power-law index b that describes the merger rate at 
large mass ratios is robustly constrained tohe b < -0.2. The 
important consequences of this finding will be discussed in 

The only regime where there is a significant difference be- 
tween results from the MS and the MS2 is where halos of 
< 100 particles in the MS are involved. In G09 we used a 
lower threshold of 40 particles, and found an upturn of the 
merger rate when the less massive halo had between 40 and 
w 100 particles. Mergers involving the same halo masses 
in the MS2 are resolved with many more particles, and no 
such upturn appears th ere. This is qualitatively understand- 
able given the finding of lWarren et a l. ( 2006) that FOF groups 
with low particle numbers are overestimated in mass. Hence, 
combining the two simulations shows that the lower thresh- 
old of 40 particles we used in G09 is too low, but that merg- 
ers involving two halos of > 100 particles are well resolved. 
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Therefore, we have excluded mergers with halos that consist 
of less than 100 particles from our global fit. This places a 
limit for the halo mass of M « 1.2 x 10" M© for the MS 
and M « 9.4 X 10^ M© for the MS2. Given that there are 
enough statistics for halos of M « 10'^ Mq in the MS and of 
M w lO''* M0 in the MS2, the largest merger mass ratios we 
can reliably probe are « 10'* with the MS and w 10^ with the 
MS2. 

2.2.2. Halo growth modes 

In the following we investigate the relation between the 
total mass growth of halos, the relative mass growth due to 
mergers and the halo merger rate. In Figure |2] the solid blue 
curves show F{< x), the fractional cumulative contribution of 
mergers to the total instantaneous growth rate of halos. Those 
contributions are summed up directly from all mergers in our 
trees. There are 6 such curves, 3 from the MS and 3 from 
the MS2, for masses spanning the range 10^ Mq to 10'** M©, 
and averaged over the redshift range 1 < z < 3 (we find only a 
weak redshift dependence ofF{< x)). Each solid curve breaks 
at some x and becomes horizontal - this is the mass ratio above 
which mergers cannot be resolved anymore, depending on the 
mass bin and simulation used^. 

The solid blue curves in Figure |2] form together a common 
envelope in the range where mergers are resolved. This en- 
velope shows that « 20% of the total growth rate comes from 
maj or mergers (1 <x < 3), mergers with 1 <x< lOcontribute 
w 30% of the mass growth, all mergers with 1 < x < 100 con- 
tribute only « 45%, and the total relative mass contribution of 
mergers even in the best resolved case, up to a mass ratio of 
10^, is no more than 60%. 

Figure |2] also shows the integration of the merger rate ob- 
tained by implementing other algorithms for merger tree con- 
struction. The blac k pluses are for th e merger rate fitting 
formula provided by IStewart et al.l (2009). Their method re- 
sults in a converging merger contribution that agrees well with 
ours'*', because they pay attention not to double-count merg- 
ers. To do that, they use a combination of the "stitching- 



oo_and splitting algorithms (see G09 and Fak houri & Mai 
l2009l for a detailed comparison of the different algorithms). 
The green symbols are for methods where some fragmen- 
tations remain in the trees. These are "snipping" (which 
is eq uivalent to not treating fragmentations at all; circles) 
from lFakhourr& Ma (200 11)" "stitching- 3" (filled circles) and 
"splitting-3" from Fakhouri & Ma (2010) (triangles), as well 
as "splitting-3" from Fakhouri et al. (2010) (filled triangles). 
The former two have an asymptotic power-law with b > 0, 
which means that the total merger mass contribution diverges 
as X increases. The latter two converge, but still show a very 
different shape from what is obtained from our trees. This 
does not mean that those methods do not conserve mass. 
Rather, as mergers with increasing x are resolved, their arti- 
ficial contribution due to fragmentation, as described in 32.11 
increases, while the compensation comes in the form of nega- 

' Here, as opposed to the case of the merger rate, we do show the contribu- 
tion of mergers with all halos, i.e. down to the resolution limit of 20 particles. 
While halos with less than 100 particles show an upturn in the merger rate, 
their influence on the mass contribution is very small, and so we include them 
in Figure |2] in order to show the full contribution of all mergers in the sim- 
ulation. The very small upturn at the higher mass ratios in some of the blue 
curves in Figure[2]are evidence for this upturn. 

'" The merger rate quantified by Stewart et all 1200% has some notable 
differences to ours in its mass ratio and redshift dependencies, but they be- 
come much less significant when the fractional cumulative contribution is 
considered, as in Figure|2] 
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Fig. 2. — The relative contribution of mergers to halo mass growth F(< x) 
as a cumulative function of mass ratio x. The soHd blue curves describe 
F{< x) for M K, lO" , 10'^, lO'"* from the MS (going up to mass ratios of 
4,400,4000, respectively) and M Si lO', lO'", lO" from the MS2 (mass ra- 
tios of 5,50,50000, respectively). The total contribution of mergers is at 
most ai 60% of the total growth rate of halos. The dashed red curves show 
the integral of the merger rate (equation (3)) using the "splitting" trees and 
normalising to the total actual growth. The green symbols show the same 
quantity, for M si 10'^ in the MS, for other m erger tree construction algo- 
rithms: "snipping" from'Fakhouri & Ma" ("2008) (circles), "stitching-3" (filled 
circles) and "splitting-3" from Fakhouri & Ma l20f3) (triangles), "splitting- 
3" fromiFakhouri et al.i >201Qi) (filled triangles), and the combined method of 
IStewart et all 120091) (vluses). 



tive contributions from smooth accretion. When mergers with 
high enough x are resolved, those methods are expected to 
give negative smooth accretion rates. 

The conclusion is that in merger trees that are built so that 
some fragmentations remain, halo mass assembly must be de- 
scribed by three components: the merger rate, the smooth ac- 
cretion rate and the fragmentation rate. Otherwise, the in- 
terpretation of "anything but mergers" as "smooth" is false. 
In ^ we use a particle-based analysis that is independent of 
the merger tree construction algorithm to show that the con- 
tribution of mergers is consistent with our fragmentation-free 
"splitting" trees. 

The finding that at most" 60% of the growth rate of halos 
is achieved via mergers with 1 < x < 10^ is already remark- 
able. But what if we had a simulation with an even larger 
dynamic range? We estimate this by using an extrapolation 
of the merger rate. The mass growth due to mergers with 
xi < X < X2 is 



/■Xi 
XI 



1 



dN„, 



'-(x,Z,M)Msmalldx, 



A'de.c-h.-uo duidx y-^i<;"^J"^smallu,^., (2) 

where M is the descendant mass and Msmaii is the mass of 
the less massive progenitor of each merger Evaluating the 

" Note that F(< x) in Figure[2]is averaged over different halos. The in- 
stantaneous value for individual halos may be very different. 
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Fig. 3 . — The ratio of the merger rate predicted from EPS to our global fit. 
This ratio is independent of redshift. The EPS merger rate is higher, which 
roughly compensates for the lack of smooth accretion in the EPS model to 
give a similar total accretion rate. 



integral of the merger rate per descendant halo requires spec- 
ifying M,„p/M, where M„,p is the main progenitor mass, since 
the merger mass ratio is defined such that Msmaii = M,np/x. 
We approximate M„p/M by -^{M^p/M), where {M„p/M) 
is the mean M„,p/M computed separately for each M in each 
of the simulations and averaged over redshift. The dashed red 
curves in Figure|2]show equation ^ evaluated between xi = 1 
and X2=x and normalised to the total actual growth. There is 
a very good agreement with the directly extracted fractions. 
This integration demonstrates that F(< x) converges at x ^ 1 
to w 60%. The convergence can be easily understood, since 
the merger rate behaves as a power-law with b = -0.3 at x ^ 1 
(see equation ([T]i). Naturally, we cannot be certain that an ex- 
trapolation is valid. Yet, for F{< x) to converge to 1, more mi- 
nor mergers are needed below the resolution limit, such that 
the asymptotic power-law index would have to be /? w -0.01 
at X > 10^. Such an index, as we demonstrate in Figure [U 
is excluded by the data at the currently available resolution 
(x < 10^), thus for the fractional mergers contribution to con- 
verge to 100%, the power-law index of the merger rate must 
change below our resolution limit. 

2.2.3. Comparison to the EPS model 

Figure [3] shows the ratio of the merger rate predicted 
from the EPS model by iLacev & Coi3 d 19931) (dashed) and 
iNeistein & Dekell (l2008l) (solid) to our global fit equation (1) 
for different masses (and independently of redshift). We iden- 
tify two regions: at jc < 100 the ratio is almost independent 
of X and ranges from « 1.6 to « 2.3 for different masses, 
while at JC > 100 the ratio is a power-law. This is because 
our fit has an asymptotic power-law index b = -Q.3, while the 
INeistein & Dekell (l2008l) merger rate has a sha llower index of 
-0. 16 ± 0.01 and that of lLacey & Coi3 (1 19931) a steeper index 
of -0.5. 

Since most of the merger-contributed accretion rate comes 
from X < 100 (Figure 0), and EPS has a w 100% higher 
merger rate in that regime, the total accretion rate of EPS 
from X < 100 equals almost the actual total accretion rate 
measured in the simulations. Since in the EPS model all 
the growth comes from mergers by construction (see however 
lAngulo & Whitd (12010)). it seems that the EPS prediction 
differs from the simulation results in two ways that roughly 



cancel each other: smooth accretion is lacking, but this is 
compensated by a boosted merger rat e. The higher merger 
rate found bv INeistein & Dekell (l2008l) in the x > 100 regime 
boosts their self-consistent EPS total accretion rate further, so 
that it overestimates the total accretion rate in the simulations 
by w 35%, with weak dependencies on mass and redshift. 

3. ANALYSIS OF PARTICLE HISTORIES 
3.1. Method 

As much as arguments exist in favour of one or the other 
algorithm for merger tree construction, some freedom is still 
left due to the complexity of the hierarchical buildup of cos- 
mic structures. As the results we presented in ^(specifically 
Figure m are algorithm-dependent, it is beneficial to perform 
an analysis that is independent of such algorithms. Compar- 
ing the results of such an analysis to the results from various 
merger trees can also serve as a tool for distinguishing be- 
tween the algorithms. In this Section we present an analysis 
of particle histories that circumvents many of the details in- 
volved in building merger trees and just relies on the identi- 
fication of structure and identification of a 'main progenitor 
trunk' for each halo. This direct particle analysis allows us to 
get a better handle on the nature of the smooth component. 

We perform this analysis on two cosmological N-body sim- 
ulations. One is the milli-Millennium Simulation that uses the 
same cosmology and has the same resolution as the MS but in- 
cludes a factor of 512 less particles in a box of 62.5/z~' Mpc on 
a side. The second is the USM Simulation (first presented in 
iMoster et al.l2010l) that uses somewhat different cosmological 
parameters that are in better agreement with current observa- 
tions (fl,,, = 0.26, r^A = 0.74, nb = 0.044, h = 0.72, n = 0.95 
and CTg = 0.77) and follows particles of mass 2 x 10^/;"' Mq 
(i.e. 4.3 times smaller than in the MS) in a box 72/i"' Mpc on 
a side. 

We distinguish between three modes of accretion: 'merger' , 
'smooth' and 'stripped'. In broad terms, we assign any par- 
ticle accreted as part of a merger event as 'merger accre- 
tion', while 'smooth accretion' is the accretion of particles 
that never belonged to a bound structure earlier than the ac- 
cretion event and 'stripped accretion' is the accretion of parti- 
cles that do not arrive as part of a halo at the time of accretion 
but were part of an identified halo at some earlier time. More 
precisely, we follow each particle p that belongs to any halo 
h at any snapshot sq to the first snapshot Sacc at which it be- 
longed to the main progenitor trunk of halo h. The halo on the 
main progenitor trunk of halo h at snapshot Sgcc is termed /!„„■• 
Note that, as we discuss below, particles may 'cycle' in and 
out of their halos, i.e. particle p does not necessarily belong 
to the main trunk of halo h at all snapshots sq > s > Sacc, but 
we are interested in the 'accretion mode' of p at the very first 
time it belonged to the main progenitor trunk of h. We then 
look for particle p in snapshot Sgcc - 1 , and tag it according to 
the following criteria. If p at Sacc - 1 belongs to a progenitor 
halo of hacc, it is tagged as 'merger accretion'. If it belongs to 
a halo that is not a progenitor of hacc, it is tagged as 'stripped 
accretion'. If p belongs to no halo at Sacc- 1, we follow it 
back through every snapshot to the initial conditions. If we 
find some earlier snapshot s < Sacc - 1 where p belonged to 
a halo, we also tag it as 'stripped accretion', otherwise it is 
tagged 'smooth accretion'. 

Further on we tag some particles as 'leaving' or 'joining' 
their halo. 'Leaving' particles are those that are not part of 
the halo's direct descendant (at snapshot s = S{)+\). 'Joining' 
particles are those that did not belong to the halo's direct main 
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progenitor (at s = sq - 1). The net growth rate of each halo is 
then given by N join- N leave- Since we also know the original 
accretion mode of each particle onto the halo, we can quantify 
the contribution of each mode to the total net growth rate. 

It is important to discuss to what extent the results from this 
analysis are expected to be independent of the merger tree on 
which the analysis is based. Let us examine the consequences 
of a fragmentation in the merger tree. The more massive halo 
of the pair, the one for which the main progenitor trunk re- 
mains intact, will be insensitive to whether the fragmenta- 
tion is cured by splitting the spuriously-connected halo or not. 
This is because a fragmentation event of a subhalo (that pre- 
viously arrived in a merger) will appear as 'leaving' particles 
tagged as 'merger mode' and therefore on average will can- 
cel out the same particles when they were 'joining'. On the 
other hand, implications do exist for the other halo of the pair, 
i.e. the smaller 'fragment', if the fragmentation is left in the 
tree. The small fragment has no progenitor and so its main 
progenitor trunk is cut and therefore the memory of its parti- 
cles as for their true origin is erased. After the fragmentation 
they will all be classified as 'stripped', because they belong 
to a halo that is not the fragment's progenitor just prior to its 
'appearance' as an independent halo. Therefore, we expect 
an artificial overestimate of 'stripped' particles if fragmenta- 
tions are left in the tree, at the expense of both other accretion 
modes. We show and discuss this effect further in Section[ 
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3.2. Results from the analysis of particles 

Figure |4] shows the fractional mass content of z = halos 
in the milli-Millennium simulation separated into the three 
modes by which each of their particles was accreted onto 
them, as indicated in the legend. The dots are for different 
halos and demonstrate the distribution around the mean val- 
ues, which are shown by curves with open circles. We see 
{green) that almost all halos are made mostly of particles that 
never belonged earlier to other halos. Plotted on top {bottom 
black) are results from the analysis of the merger trees, which 
are obtained by summing up all the mass accreted via mergers 
along a halo's main progenitor trunk. An important point to 
take from Figure |4] is the agreement between the two analy- 
sis methods, shown by the red and bottom black curves. The 
small difference between the two is expected because some 
particles that arrive via mergers later leave the halo and do 
not contribute to its mass at z = (see also below). 

The trend seen in Figure|4]is clearly an effect of resolution: 
halos closer to the resolution limit have gained most of their 
mass by unresolved accretion, which cannot be distinguished 
between the different components and this mass is identified 
as 'smooth' . Nevertheless, the beginning of a saturation of the 
mass fraction arriving from mergers is apparent as halos are 
better resolved. We suggest that the 'true' value for all halos 
is the saturated value that is seen for high-mass halos. Addi- 
tional evidence for the trend being an effect of resolution is the 
fact that for the same halo mass, the MS2 halos, which are bet- 
ter resolved, show a weaker trend {top black) than the results 
from the MS {bottom black). We do not have enough statis- 
tics to quantitatively constrain the saturation of this quantity 
due to low halo numbers at high masses, but these results are 
consistent with our results on the accretion rate ( see below), 
where the saturation is statistically robust 02. 2.21) . 

In Figure |5] we show the main result of the particle analy- 
sis, the fraction of the net instantaneous accretion that is as- 
sociated with the three different modes, as a function of halo 
mass and redshift. We observe no dependence on redshift. 
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Fig. 4. — The fraction of the particles in z = halos originating from the 
three accretion modes: 'merger' (red), 'stripped' (blue) and 'smooth' (green), 
versus halo mass. The dots are for different halos and demonstrate the distri- 
bution around the mean values, which are shown by curves with open circles. 
Overplotted are corresponding curves from the analysis of the merger trees of 
the MS (lower black curx'e) and MS2 (upper black curve), which indicate the 
total mass accreted via mergers onto the main progenitor trunk normalised to 
the z = halo mass. The error bars include 68% of the halos around the mean 
(see McBride et al. 2009 for a more detailed description of the distribution 
around the mean for major mergers). The MS curve from the merger trees 
is similar, yet slightly higher, than the 'merger' component from the particle 
analysis. This is expected since some particles that arrive via mergers leave 
(smoothly) and never come back, which is information that is not included 
in the merger trees alone. The MS2 curve is significantly higher where MS 
halos ai'e barely resolved and MS2 halos ai'e well resolved, but the difference 
between the two curves decreases towards higher halo masses, as expected. 
Note that the mass content of halos is only slightly affected by the proxim- 
ity to the end of the simulation and the inability to split some subhalos out. 
This is because the mass ratio regime that is mostly affected by this effect is 
subdominant in mass contribution. 



even down to z = 0, despite the inflated merger rate that is 
due to the proximity to the end of the simulation. This means 
that even if some of the increase of the minor merger rate at 
z < 0.5 is real, it is not high enough to make a significant con- 
tribution to the total mass accretion (see also Figure|9]below). 
The mass dependence is again a resolution effect. As halos are 
better resolved, the mergers they undergo are better resolved 
and so the contribution of the 'merger' component increases. 
Although we again cannot make a robust fit and extrapola- 
tion from these results, we find it reassuring that the particle 
analysis gives a consistent result with that of the "splitting" 
merger trees {dashed black curves, which cover the range of 
the dashed curves in Figure|2]i. 

It is worth noting that the curves in Figure|4]are shifted from 
those in Figure|5]by ?» 0.5 - Idex toward higher masses. This 
is expected, as halos accrete most of their final mass M^-q, by 
definition, when M^ > 0. IM—o. This means that the fractions 
of the different accretion modes in the final halo mass M— o 
(FigurelU correspond to their fractions in the accretion at M > 
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Fig. 5. — The fraction of net accretion (i.e. 'joining' particles minus 'leav- 
ing' particles), as a function of halo mass, that belongs to each of the three 
modes 'merger' (red), 'stripped' {blue) and 'smooth' (green). The results are 
shown for three different redshifts and are independent of it. As higher mass 
ratios are resolved, more of the accretion is in the 'merger' and 'stripped' 
modes, but the merger mode contribution is consistent with the saturation in- 
ferred from the merger trees (shown by the dashed black curves that cover the 
range of the dashed curves in Figure[2). This indicates that the contribution 
of mergers may not increase further. 



0.1M.=o (Figure©. 

3.3. The cycle of particles in and out of halos 

From Figures |4] and |5] we can also learn that the fraction 
of the accretion in the 'stripped' mode is subdominant to that 
in the 'smooth' mode. That is, what we could only interpret 
as 'non-mergers' from the analysis of the merger trees, can 
now be shown to consist of particles that never belonged to 
another halo prior to their accretion. In fact, the 'stripped' 
mass is consistently « 1/3 of the mass in the 'merger' mode. 
Conservatively, the 'stripped' component is our uncertainty, 
because our analysis does not indicate how long it has been 
stripped and what part of it comes from the vicinity of ap- 
proaching subhalos. Indeed, some of these 'stripped' par- 
ticles arrive as mergers into low-mass halos and after being 
stripped from them, are re-accreted into more massive halos 
and tagged there as 'stripped' . Thus, some of the mass appear- 
ing in the merger trees as 'merger mode' is transferred into 
'stripped mode' of more massive halos in the particle analy- 
sis. This is the reason the red curves are somewhat lower than 
the black dashed ones. 

We can learn more about this 'cycle' of particles through 
different halos from Figure |6l There it is shown that the rate 
at which particles join and leave their halo is similar to, or 
even higher than, the net growth rate, for each of the differ- 
ent modes. As a quantitative example, shown in Figure|6] the 
rate of 'smooth' particles joining their halos at z w is « 4 
times higher than the net growth rate due to smooth accre- 
tion, i.e. only « 33% higher than the 'leaving' rate of parti- 
cles that previously arrived smoothly. These numbers drop 
toward higher redshift, where the 'cycle' is less significant, 
e.g. at z « 2 the values are lower roughly by a factor of 2 
compared to those shown in Figure |6l Note that all 'leaving' 
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Fig. 6. — The z = cycle of 'joining' and 'leaving' particles, shown by the 
ratio of the gross inflow rate ('joining' particles) to the net growth rate ('join- 
ing' minus 'leaving'). Both the gross and net rates are calculated separately 
for each accretion mode, which are shown by the different curves in indicated 
in the legend. The most striking result here is that both the gross mass gain 
and the gross mass loss are comparable to, or even larger than, the net halo 
growth rate. 



particles leave their halos smoothly, i.e. not as part of a bound 
subhalo, as such events have already been cleaned at the time 
of merger tree construction. Thus the red curve in Figure |6] 
shows that there are many particles that arrive via mergers 
and then stripped off of their subhalos inside the main halos 
and later leave the main halo 'smoothly'. It seems reasonable 
to suggest that this cycle is driven, at least partly, by fluctua- 
tions of particles that reside close to the halo boundary in and 
out of the region defined as the halo by the FOF algorithm. 
This is the reason we focus throughout the paper on the net 
growth rate. We leave a more detailed study of this cycle to 
future work. 

3.4. Alternative trees and halo definitions 

In Figure |7] we compare the merger contribution to 
the accretion rate from our "splitting" and "snipping" 
(Fakhouri & Ma 2008) trees. The "snipping" algorithm is the 
most extreme case of leaving all fragmentations in the trees so 
that in the merger tree analysis the merger contribution does 
not converge {pluses, same as in Figure |2]l. In contrast, the 
"snipping" merger contribution does not exceed 60% when 
the particle histories are used (dashed). In fact, it is somewhat 
lower than in th e "sp litting" case (solid, asterisks) because, as 
we describe in ^3.11 the role of 'stripped' particles is overes- 
timated at the expense of the other accretion modes. This is 
because the particles of 'snipped' halos, which lose their main 
progenitor trunks, are tagged 'stripped' instead of their orig- 
inal accretion modes. To verify this we divide the 'stripped' 
particles into two categories: 'just stripped' and 'stripped in 
the past' . The former category includes particles that at Sg^c - 1 
belong to a halo that is not a progenitor of h„cc, and the latter 
category includes particles that were found to belong to a halo 
at an earlier time in the past. We find that 'just stripped' parti- 
cles are negligible in the "splitting" trees but they become very 
significant, especially for low mass halos, in the "snipping" 
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Fig. 7. — Comparison of the contribution of the different accretion modes 
at z = 0.3 in the "splitting" and "snipping" trees from the USM simulation. 
In the "snipping" trees (dashed) the contribution of mergers is not higher 
than in the "splitting" trees (solid, asterisks), unHke the results of the merger 
tree analysis (pluses, same as in Figure |2j. This demonstrates that merger 
trees with fragmentation events are not appropriate for studying halo growth 
modes, as they are not self-consistent. In fact, in the particle analysis the 
"snipping" merger contribution is somewhat lower than the "splitting" merger 
contribution due to an artificial effect caused by the snipping of the main 
progenitor trunk, see details in the text. For the "splitting" trees, there is 
very good agreement between the USM simulation (solid, asterisks) and the 
mini-Millennium simulation (Figure |5). 



trees, because the particles of a 'snipped' halo liacc are identi- 
fied at Sacc - 1 in the artificially-connected FOF group, which 
is not a progenitor of hacc Hence, we see that the results of 
the analysis we present in this Section are not entirely merger 
tree-independent. The results for trees with remaining frag- 
mentations are different depending on the analysis method 
(merger tree only or particle history analysis) and both suffer 
from artificial effects because the trees are not self-consistent. 
On the other hand, the analyses based on the "splitting" trees 
are consistent with each other (as already shown in Figure|5]). 
Finally we explore the sensitivity of our results to the halo 
definition by repeating our analysis with a different defini- 
tion for a halo. The first step of structure identification, the 
FOF group finder, is run with different linking lengths of 
b = 0.25 (b = 0.15), corresponding approximately to struc- 
tures with overdensities of 100 (500) rather than 200. This re- 
sults in FOF groups that are on average w 15% more (w 25% 
less) massive. We then run SUBFIND and our splitting al- 
gorithm as before. The halos in our final catalogues are also 
more (less) massive compared with our original catalogues 
by similar factors. This is because SUBFIND is run on each 
FOF group separately and is restricted to working on particles 
within FOF groups only '^. Figure|8]shows the contribution of 
the different accretion modes to the halo mass at z « 0.3 from 
the USM simulation for the three different FOF group defini- 
tions. The results agree very well with each other and with 

'^ In the b = 0.25 case, where particles have been added to the FOF groups 
compared with the b = 0.2 case, they are typically infalling onto their halos, 
i.e. they have negative energies relative to the halo centers and are therefore 
considered bound by SUBFIND and included in the new halo catalogues. 
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Fig. 8. — The contribution of the different halo growth modes to the z = 0.3 
halo mass from the USM simulation. The results with the standard FOF group 
finder (solid) agree very well with the results from the milli-Millennium 
simulation (Figure |4). The results based on different halo definitions (FOF 
groups with b = 0.25, dashed and with b = 0.15, dotted) are also virtually 
unchanged. 



the results presented above from the Millennium Simulations. 
This agreement reflects the fact that the change of halo defini- 
tionis roughly mass independent (Jenkins et al. 2001; Whitg 
I2OO2I) . and so both the merger mass contribution and the to- 
tal halo mass change in a similar way. It then follows that 
the smooth accretion mass contribution increases by the same 
factor and the relative contributions do not change. Note, 
however, that using a mass definition that changes the mass 
of halos as a function of mass (e.g. Cuesta et al. 2008) could 
lead to different results. For example, Warren et al.. d2006b 
suggested a correction to FOF group masses that reduces the 
mass of poorly resolved (low mass) FOF groups, an effect 
that would work in the direction of an even higher contribu- 
tion from smooth accretion. 

There exist a few alternatives to the FOF algorithm for iden- 
tifying dark matter halos. In this paper we use only the FOF 
algorithm (with different linking lengths, as described above), 
but we believe that other halo definitions will not change the 
m ain results. One su pporting evidence for that is the work 
of IStewart et al.l ( 2008,). who used th e 'Bound Density Max- 
ima' algorithm (iKlypin et al.l Il999h for halo identification, 
and found results that agree well with ours. Another alterna- 
tive would be to use 'Spherical Overdensity' halos, i.e. defin- 
ing the spherical region inside the virial radius as the halo. 
Unless FOF groups are organised in a way that the smoothly 
accreted mass is outside the virial radius and the merger ac- 
creted mass is inside, significant deviations from our results 
should not arise. The case is probably the opposite, because 
when FOF groups have very aspherical shapes it is mainly be- 
cause a few substructures are connected together, i.e. it is not 
that the outskirts of FOF groups constitute mainly of smoothly 
accreted material. 



4. IMPLICATIONS FOR GALAXY FORMATION 

Many theoretical models focus on the role of major mergers 
in galaxy evolution. In this paper we have shown that it is ac- 
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tually minor mergers {x> 10) and smooth accretion that dom- 
inate halo growth by accounting for w 70% of the accretion 
rate onto halos. Therefore, halos grow mainly continuously 
rather than in bursts with short duty cycles. This mode of ac- 
cretion is more favourable for disk formation at all redshifts, 
and may in particular help in understanding galaxies at high 
redshift that show extr eme star-formation rates but no signs 
of major mergers (e.g. Genzel et al.ll2008l; iGenel et al.ll2008l; 
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If there is a truly smooth accretion component, simu- 
lations from cosmological initial conditions naturally in- 
clude it. Indeed, our result is consistent with hydrody- 
namical simulations that show that most baryonic accre- 
tion onto galaxies does not arrive in the form of mergers 
jMurali et al." '2002"; 'Sem eUn & CombesI 120051: iMaller et all 
I2OO6; Dekel et al. 2009a). Our findings suggest that as the 
resolution of future simulations increases, the mass contri- 
bution of small halos to the formation of galaxy-size halos 
will hardly increase further. This implies that « 40% is a 
strong lower-limit on the mean fraction of pristine IGM gas 
in the baryons accreted onto halos of any given mass or red- 
shift, since this smooth gas was never bound to any subha- 
los and is therefore not expected to have formed stars or to 
have become significantly enriched with metals, regardless of 
the star-formation efficiency and history of the baryons in the 
merging halos. It may also be expected that this w 40% did 
not experience feedback from star-formation in smaller halos, 
and is likely "cold" when finally accreted onto a halo, with 
T > lO'^K set by photoheating of IGM gas. 

In the context of the "cold flow" mode of gas accretion onto 
galaxies (Birnboi m & Dekel 2003; Dekel & Birnboim 200^ 
iKeres et aLll2005l) . our results suggest that roughly half of the 
incoming gas is not clumpy, thus processes like tidal stripping 
and dynamical friction are irrelevant for it'^. Theoretical or 
semi-analytical galaxy formation models, for example EPS- 
based models, may benefit from taking the smooth accretion 
component into account, as many properties of galaxies, such 
as size and morphology, depend on the way they accrete their 
baryons. 

Even if an extrapolation of the power-law index of the 
merger rate beyond x= 10'' is not valid and mergers of higher 
mass ratios do make up for the "missing" « 40%, our results 
are still significant. Since the MS2 resolves all halos with 
r,„> > lO'^K at z < 3 (M > 1.9 X 10^ M©), we do resolve most 
of the accretion of halos that have experienced star- formation, 
that is if star- formation is prohibited in halos with ri,,> < 10"* A " 
(e.g. iR ees 1986: Efstathioul ll992l: lOkamoto & Freiifl 120091: 
iHoeft & Gottloeber.2010.) . Current models of galaxy forma- 
tion derive, assume or require strong suppression of cool- 
ing and/or star-formation below a virial temperature thresh- 
old that is even higher than lO'^K (e.g. iBouche et al.l 120091 : 
lKravtsovl l2010, and references therein). The fraction of the 
mass that arrives either smoothly or by accretion of small ha- 
los can be read off Figure |2] as the complement of the merger 
contribution. However, for the convenience of the reader we 
explicitly show in Figure |9] the fraction of accretion that ar- 
rives smoothly, for different values of a threshold mass be- 
low which accretion of baryons is assumed smooth. The re- 
sults shown in Figure |9] do not depend on any extrapolation 
but are based directly on the MS2 (since all M,/, we use are 

'^ Although it is possible that the gas becomes clumpy inside the halo on its 
way t o the galaxy due to different instabilities (e.g. Field 1965; Burkert & Lin 
l2000tlMaller & Bullockii2004ciKeres & Hemquislii2009i : iBimboina2009ii 
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Fig. 9. — The fraction of mass accretion that arrives "more smoothly than" 
M,i,, i.e. by accretion of halos with M < M,/, plus smooth accretion. This 
fraction is shown as a function of halo mass, and for different values of M,;,. 
These results are from the MS2 and are based directly on the merger trees, 
without using the fitting formula Equation [T]or any extrapolation below the 
simulation's resolution limit. The results are plotted for accretion at j f» 1 
(top) and z ~ (bottom). The resemblance of the z ~ results to those at 
z ft! 1 highlights the fact that despite of an overestimation of the minor merger 
rate at z < 0.5 in the "splitting" trees, the contribution of minor mergers to 
the total mass growth is still very small. 



resolved by the MS2). For example, by assuming that pre- 
heating from the cosmic UV background evaporates halos be- 
low M,h « 10^ M© (green), we infer that 10" M© halos at 
z « 1 (z « 0) accrete w 55% (w 65%) of their baryonic mass 
in smooth T > IQ'^K gas. Since 10" M© halos do not form 
a stable virial shock and the cooling times are short, over- 
efficient galaxy formation in such halos is probably prevented 
by strong baryonic feedback processes. 

5. COMPARISON TO PREVIOUS WORK 

Merger trees are regularly treated in the literature such that 
each node (whether a halo or a subhalo) never has more than 
one descendant. This constraint is motivated by the idea of 
hierarchical buildup of structure, but in practice, when stan- 
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dard halo definitions are used, fragmentation events are com- 
mon too. To treat the formation of dark matter halos self- 
consistently, there are two alternatives. First is to allow for 
more than one decendant and consequently quantify both the 
merger rate and the fragmentation rate. Alternatively, the 
trees may be rearranged so that fragmentations do not ex- 
ist. In this paper we use a variant of the second possibility, 
namely our "splitting" algorith m'"^. Our result agre es with 
that of Ste wart et al. (2008) and I Stewart et al J (120091) . as de- 
scribed in ^2.2.21 since they use fragmenta tion-free merger 
tree s, but d isagrees with the results of Fakhouri & Mai (l2008h 
and lFakhou ri & Ma (2009^, who d o not. 

Recently Fakhouri &Mal (1201 Ol) found that the relative im- 
portance of mergers to the growth of halos in the MS cor- 
relates with large-scale environment. This, as they suggest, 
is further compelling evidence for the true diffuse nature of 
non-merger halo growth, even if the quantitative proportions 
t hey find from their tree s differ from ours. 

iGuo & White! (l2008h investigated the relative growth of 
FOF groups in the MS via major mergers, minor mergers 
and 'diffused particles'. Their main findings that are differ- 
ent from our results are: (i) at large enough masses (M > 
3 X 10'^ Mq) merger accretion dominates over smooth accre- 
tion, (ii) the relative role of merger accretion, for a given 
mass, is larger at lower redshifts. Both differences occur 
because Guo & White (2008) consider "unprocessed" FOF 
groups, i.e. implicitly use 'snipping' trees. They consider only 
instantaneous mass gain due to mergers but not instantaneous 
mass loss due to fragmentations, which becomes more impor- 
tant at low redshift and for high mass (better resolved) halos. 

Angulo & White (2010) studied halo formation histories in 
the framework of EPS assuming a finite dark matter particle 
mass. They found that if the dark matter particle is assumed 
to be a 100 GeV neutralino then Milky Way type halos are 
expected to have w 10% smooth accretion, which is due to 
a minimum halo mass imposed by the free streaming of the 
dark matter particles. They compare their results to an N- 
body simulation, and find that up to a mass ratio x = 500 the 
contribution of mergers is k, 45%, in good agreement with 
our results. For the same mass ratio, they find the spherical 
collapse EPS model to give « 70% accretion from mergers, 
in rough agreement with our results as well. The significant 
contribution they find for mergers with x > 500 (« 25% in 
the spherical collapse model) is due to the higher power-law 
index of the merger rate in EPS, as shown here in Figure [3] 
If instead the power-law index appropriate for very large x is 
the one we find from the resolved regime in the Millennium 
Simulations, the contribution of mergers to halo growth prac- 
ti cally vanishes even b efore reaching the free streaming mass. 

iMadau et all (12008 ') report that less than 3% of the final 
mass of the Via Lactea halo, which is a Milky Way-type halo 
simulated at high resolution, was accreted smoothly. This ap- 
pears to stand in sharp contrast with our results, and there are 
a few possible explanations for this disagreement. 

• We do find Milky Way-type halos that have < 3% of 
smooth accretion, but they are only approximately one 
in a thousand halos. It could then be that the Via Lactea 
halo is an untypical halo (« 3a) in this respect. This is 
a possible but undesired 'last resort' explanation, which 
we do not need to invoke, given the more likely expla- 

''* In G09 we found that there is a good physical basis for our specific 
choice of a "splitting" algorithm, since fragmentations usually occur before 
the halos have had significant dynamical interaction. 



nation below. 

• It could be that the higher resolution available for the 
Via Lactea halo is responsible for the difference. We 
suggest that this is not the dominant reason for the dif- 
ference based on the following argu ment. It can be read 
from Figure 4 in lMadau et al.l (2008) that halos of peak 
velocity ratio < 45 have contributed to the Via Lactea 
halo 97% of its final mass, which is all of the accretion 
that is associated with mergers. Such velocity ratios 
correspond approximately to mass ratios x < 45^ sa 10^, 
which can also be resolved for the most massive ha- 
los in the MS2. In other words, it seems that these 
are not the higher resolution mergers in the Via Lactea 
halo that contribute the "additional" (when compared 
with our results) w 37% mass, but rather mergers that 
we are able to resolve. Moreover, we can consider just 
the total fraction reported by Mada u et al.l (12008). The 
maximal mass ratio resolved by the Via Lactea simula- 
tion is jc sa 5 X 10^, and our best resolved mass ratio is 
jc « 10^. For this range 10^ < x < 5 x 10^ to account for 
« 97% -60% = 37% of the mass accretion, the power- 
law index of the merger rate would have to be /? w 1 .5 at 
10^ < X < 5 X 10''. Compared with the value b « -0.3 
we find atx < 10^, this is an unlikely abrupt change. 

• We beheve the most probable explanation is that differ- 
ences i n analysis methods cause the results to be so dif- 
ferent. iMadau et al.l (12008 ) sum up all the peak masses 
(over their full formation histories) of all halos that have 
merged into the main progenitor trunk and compare it 
to the mass inside r2oo at z = 0. They do not account for 
the fact that some fraction of the mass of the merged 
subhalos ends up outside of r2oo at z = 0, and so they 
overestimate the mass contribution of mergers. 

To understand this issue with greater certainty, other high res- 
olution simulations of individual halos should be further ex- 
amined. To this end, J. Wang et al. 2010 , in preparation, anal- 
yse the high resolution Aquarius halos dSpringel et al.ll2008l) 
with a method close to the one we use in ^ and find re- 
sults that are more consistent with ours than with those of 
IMadau etal.1 (120081) . 

6. SUMMARY 

In this paper we calculate the merger rates of dark matter 
halos and we investigate the role of smooth accretion versus 
mergers to their growth. We extract the merger rates and ac- 
cretion histories from the Millennium and Millenium-II simu- 
lations, combined with two additional, smaller, cosmological 
simulations. We use the "splitting" merger tree construction 
algorithm described in ^2. H and in G09, and verify its reliabil- 
ity by reproducing our results by following individual particle 
histories alone, independent of merger tree construction algo- 
rithms. 

We find that the contributions of all resolved mergers (up 
to mass ratios « 10'') to the total growth rate of halos do not 
exceed 60%, regardless of halo mass and redshift. Most of 
the merger contribution comes from small mass ratio ("ma- 
jor") mergers (e.g. 1 < x < 10 contribute « 30% of the to- 
tal growth), while "very minor" mergers add very little mass 
(e.g. 10^ < X < 10^ contribute just a few percent to the total 
halo growth). We find that the power-law index of the merger 
rate is such that if the merger rate is extrapolated beyond the 
maximum resolved mass ratio « 10^, the total contribution of 
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all mergers saturates at « 60%. This suggests that a signif- 
icant mass fraction of halos may be accreted in a genuinely 
smooth way. 

Our results have important implications for galaxy forma- 
tion models and the modes of baryonic accretion. If w 40% of 
the dark matter that is accreted onto a halo was never previ- 
ously bound in any merging smaller halos, then at least « 40% 
of the baryons must also be accreted smoothly - as gas that 
was never previously heated by feedback processes or con- 
verted to stars. For the baryons, this « 40% is a strong lower 
limit since halos with ry,> < 10'^ K likely cannot retain their 
gas. The common assumption that halos with T^r < IQ'^K 
cannot retain their gas also makes our results insensitive to 
the limited resolution of the simulations we use, because we 
resolve all halos above this limit (at z < 3). The implication is 
that a very large fraction of the baryonic matter falling into a 
halo must be pristine "cold" IGM gas, with T > 10*K set by 
IGM photoheating. This gas is not expected to have formed 
stars or to have become significantly enriched with metals, no 
matter what the star-formation efficiency and history of the 



baryons in any of the merging halos. 
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